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Abstract 



A three-dimensional approach based on momentum vectors as variables for 
solving the three nucleon Faddeev equation in first order is presented. The 
nucleon-deuteron break-up amplitude is evaluated in leading order in the NN 
T-matrix, which is also generated directly in three dimensions avoiding a 
summation of partial wave contributions. A comparison of semi-exclusive 
observables in the d(p, n)pp reaction calculated in this scheme with those 
generated by a traditional partial wave expansion shows perfect agreement at 
lower energies. At about 200 MeV nucleon laboratory energies deviations in 
the peak of the cross section appear, which may indicate that special care is 
required in a partial wave approach for energies at and higher than 200 MeV. 
The role of higher order rescattering processes beyond the leading order in 
the NN T-matrix is investigated with the result, that at 200 MeV rescattering 
still provides important contributions to the cross section and certain spin 
observables. The influence of a relativistic treatment of the kinematics is 
investigated. It is found that relativistic effects become important at projectile 
energies higher than 200 MeV. 
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I. INTRODUCTION 



During the last two decades calculations of nucleon-deuteron (Nd) scattering based on 
momentum space Faddeev equations [1] experienced enormous improvement and refinement. 
It is fair to state that below about 200 MeV projectile energy the momentum space Faddeev 
equations for three-nucleon (3N) scattering now can be solved with very high accuracy for 
the most modern two and three nucleon forces. A summary of these achievements is given 
in Ref. [2]. The approach to 3N scattering described in Ref. [2] is based on using angular 
momentum eigenstates for the two- and three-body system. For low projectile energies 
this procedure is certainly physically justified due to arguments related to the centrifugal 
barrier. However, to probe the strong interaction at shorter distances one has to go to higher 
projectile energies, where the algebraic and algorithmic work carried out in traditional partial 
wave (PW) decomposition can be quite involved. A more crucial hurdle is posed by the fact 
that in three-nucleon scattering (3N) calculations for projectile energies of a few hundred 
MeV the number of partial waves needed to achieve numerical convergence proliferates, 
and limitations with respect to computational feasibility and accuracy are being reached. 
It appears therefore natural to abandon PW representations completely and work directly 
with vector variables, if one wants to calculate 3N scattering at higher energies. As an 
aside, the use of vector variables is common practice in bound state calculations based on 
variational [3] and Green's Function Monte Carlo (GFMC) methods [4], which are carried 
out in coordinate space. 

Momentum space calculations within the Faddeev scheme which did not employ a P W 
decomposition were first carried out for a system of three bosons [5,6]. Here the momentum 
space Faddeev equations were solved for the bound as well as the scattering state. In 
this work we want to employ realistic nucleon-nucleon (NN) interactions in a 3N scattering 
calculation. This means we have to incorporate spin degrees of freedom into the formulation 
of the Faddeev equations. Since the input to any Faddeev calculation is the solution of 
the Lippmann-Schwinger (LS) equation for the two-nucleon T-matrix, we start from the 
formulation of NN scattering developed in Ref. [7]. There we chose an approach based on 
the total helicity of the NN system as spin variable. From our point of view this is the 
preferred starting point to later progress to the 3N system. 

In this work we consider the first term of the multiple scattering series built up by the 
Faddeev equations, rather than solve the full Faddeev equations for three nucleons, and con- 
centrate on semi-exclusive break-up observables. Of particular interest are the polarization- 
transfer coefficients in the (p,n) charge exchange reaction on the deuteron, which recently 
have been measured at IUCF with a projectile energy of 197 MeV [8] and at RCNP with a 
projectile energy of 346 MeV [9]. Since these measurements are carried out at 'intermediate 
energies', one can speculate that it may be sufficient to consider only the first order term 
in the multiple scattering series. Furthermore, since the projectile energies are high, we will 
consider relativistic effects as far as the kinematics is concerned. 

In Section II we formulate the Nd break-up process in a three-dimensional (3D), nonrel- 
ativistic Faddeev scheme. We derive the leading term of the full Nd break-up amplitude, 
where NN T-matrix elements are given in the momentum- helicity basis defined in Ref. [7]. 
In Section III we introduce relativistic kinematics into this formulation. We will not consider 
a boost of the NN T-matrix [10] nor Wigner's rotations [11] of the spin. The observables for 
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the (p,n) charge exchange reaction are introduced in Section IV. In Section V we present 
and discuss our calculations for the (p,n) charge exchange reaction in the proton-deuteron 
(pd) break-up process. Here only the outgoing neutron is detected after the break-up. We 
present calculations of spin averaged differential cross sections, neutron polarizations, proton 
analyzing powers, and polarization-transfer coefficients at different energies. We also com- 
pare our calculations with traditional PW calculations. Finally, we summarize in Section 
VI. 



II. THE NONRELATIVISTIC ND BREAK-UP AMPLITUDE 

In the Faddeev scheme the operator Uq u11 for the Nd break-up process is given as [2] 



u£ ul1 = (1 + P)T F . 



(2.1) 



Here T F is the Faddeev operator obeying the Faddeev equation [1] for the break-up process 
of three identical particles, 

T F = TP + TG PT F . (2.2) 
The operator T stands for the NN t-matrix, and P is a permutation operator defined as 

P = -P12-P23 + Pt&Pts- (2-3) 

The free three-nucleon (3N) propagator is given by Go- The matrix elements of the break-up 
amplitude £7 /u "(p,q) of Eq. (2.1) are defined as 



U$ ull (p, q) = ( pqm sl m S 2m s3 T 1 r 2 T3 



U, 



full 



= < pqm sl m s2 m s3 T 1 T2T 3 



(1 + P)T F 



q m sl r 1 v d 



(2.4) 



where 



\pqm s im s2 m sZ T l T2T 3 ) = \qm sl T 1 )\pm s2 m s3 T 2 T 3 ) 

is the final, not-antisymmetrized free 3N state. The quantities m si ,Ti ('< 
spins and isospins of the three nucleons. The initial state, in which only the deuteron state 
^Y d ") i s antisymmetrized, is given by 



(2.5) 

1, 2, 3) are the 



n m° T°\I/ Md \ — 



(2.6) 



The index indicates the projection of total angular momentum of the deuteron along 
an arbitrary z-axis, rnP sl ,T® are the spin and isospin of nucleon '1' acting as the projectile. 
Without loss of generality nucleon '1' is singled out as projectile, while the other two nucle- 
ons, '2' and '3' form the two-nucleon (2N) subsystem, i.e. the deuteron in the initial state. 
Jacobi momenta p and q are used to describe the 3N kinematics in the final state, 
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q 



2 ( k 2 - k 3) 



ki - - (k 2 + kg) 



ki - -k te6 . 



(2.7) 
(2.8) 



Here the kj's (i — 1, 2, 3) represent the laboratory momenta of the three nucleons. Defining 
ki ab as the laboratory momentum of the projectile and applying momentum conservation 



Hab 



k x + k 2 + k 3 , 



(2.9) 



leads to Eq. (2.8). In the initial state q is the relative momentum of the projectile to the 
target deuteron and is related to kz a & as 



Qo — -k; a fc- 



(2.10) 



For clarity of description we will in the following denote the break-up amplitude as Uq u11 (p, q) 
and suppress all other quantum numbers. 

In this work we only want to consider the leading term of the full break-up operator 
Ui ul1 . This means we only consider the leading term in the Faddeev operator T F of Eq. (2.2) 
and define the break-up operator Uq in first order in the NN T-matrix as 



Uq = (1 + P)TP. 



(2.11) 



The matrix elements of £/ (p, q) with respect to the final and initial states from Eqs. (2.5) 
and (2.6) are then given as 



^o(p,q) = \ P<irn sl rri S 2rn s3 T 1 T2T 3 
= (pqm sl m s2 m s3 T 1 r 2 T 3 



U 



(1 + P)TP 



(2.12) 



From now one we mean by the Nd break-up amplitude (or break-up amplitude) the matrix 
element U (p, q) given in Eq. (2.12) and by the full Nd break-up amplitude (or full break-up 
amplitude) the matrix element Uq u11 (p, q) given in Eq. (2.4). 

The break-up amplitude C/ (p, q) from Eq. (2.12) is composed out of three terms, 



U (p, q) = t/ (1) (P, q) + £/<T(P, q) + Uf>(p, q), 



(2)/ 



r(3)/ 



with 



Uq\p,(1)= (pqm sl m s2 m s3 T 1 T2T 3 

Uo 2) (P,q) = _ ( Pq^sl^ S 2^ S 3T-lT"2T"3 

Uo 3 \p,q) = ^ { pq/»,i/".,2"'.s:;-j-2 7-:i 



TP qomW^" 
P 12 P 23 Tp\q m° sl T 1 ^ d 
P 13 P 23 TP q m sl r^ d 



(2.13) 

(2.14) 
(2.15) 
(2.16) 



Here the final free 3N state are labeled '1', meaning that nucleons '2' and '3' form the 2N- 
subsystem. The spin and isospin quantum numbers must be read in the order 123. Applying 
the permutation operator P\ 2 P 23 to the final state of £7o^(p,q) given in Eq. (2.15) leads to 
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t/ (2) (p,q) 



P23Pi2Pqrn sl m s2 m s3 TiT2T 3 



TP 



pqm s im s2 m S 3rir2r 3 



TP 



(2.17) 



This is now the final state where nucleons '1' and '2' form the 2N subsystem. Accordingly, 
the spin and isospin quantum numbers associated with the three nucleons read in the order 
312. In order to have the same final states as Uq 1 \p, q) in Eq. (2.14), we need to transform 
the final state such that nucleons '2' and '3' form the 2N subsystem. This transformation is 
achieved by the following relation among the Jacobi momenta [12] 



1 3 
Pi = -^Ps - ^3 



1 

qi = Ps - ^qs, 



where the labels 1 and 3 indicate the nucleon being singled out. This leads to 



P - 7q P - 7:q m s2m s3 m sl T 2 T 3 Ti 



Using another relation between Jacobi momenta [12] 



TP 



(2.18) 



(2-19) 



Pi = -2P2 + ^2 



qi = -P2 - 7^2' 



(2.20) 



we can obtain uff* (p, q) in a similar fashion as 

^"o 3) (P> q) = i ( (~^P + ^q) (-P - ^q) m s3 m sX m s2 T 3 T X r 2 



TP 



q m sl r 1 w d 



(2.21) 



Since Uq 2 \p, q) and uff\p, q) differ from UQ 1 \p,q) on ly i n their variables, it is sufficient 
to work out an expression only for one of them, which we choose to be Uq 1 \p, q). For 
calculating subsequently U^f 1 (p, q) and ujf 1 (p, q) one only has to perform the following 
replacements, 



for £/d 2) (p,q) : {T, ™}{1,2,3} -> {^^}{2,3,l} 



"2 P "4 q ^^P-2 q 
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for f/ (3) (p,q) : {r,m} {li2 ,3} ^ {r,m} {3ili2 } p -> --p + -q q -> -p - -q. 



(2.22) 
(2.23) 



For calculating £/o^(p, q) we start by inserting the following completeness relation for 
the free 3N system 



^ I dp I dq\pqm sl m s2 m s3 T 1 T 2 T 3 }(pqm sl m s2 m s3 T 1 T 2 T 3 \ 

m „ >-» m „ o v v 



(2.24) 



m sl m s2 m s3 
T 1 T 2 T 3 



twice into Eq. (2.14), which leads to 

^o (1) (P,q) = Yj d P/ ^q / (pq^i^2m S 3Tir 2 T 3 |T|p / q / m;i^2K3 T i r 2 T 3> 



si si s>5 

t't't' 
1 2 3 
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x £ / dp" f rfq" (p^m' sl m' s2 m' s3 r[r^\ P |pV<m>M^> 

II III!" ** 

m m m. 



m , m _ m ' 
si si sA 
^11 ll II 
T l T 2 T 3 



,0 _0 lT ,M d \ 



x (p q m^m^m^r! r 2 r 3 qom s i^i^d 7 

5] / dp' (pm S 2m s3 r 2 r3 |T(£ p )| p'm'^m'^r^) 

J r 2 T 3 

E / V (p'm^m'^| (qm sl n |P| qom^ ) IpXV^') 
„n ^.11^.11 J 



m 's2 m 's3 T 2 T 3 



m 's2 m 'U T 2 r 'i 

x (p"m>"3^' 



(2.25) 



In arriving at the last equality we used the fact that T acts only in the two-particle 
subsystem together with momentum space properties of the initial state. 

The NN T-matrix is calculated at a center of mass (cm.) energy E p of the 23-subsystem 



E„ 



P ~ = f-(?o 2 - Q 2 ) + E d , 
m Am 



[2.2Q) 



which does not necessarily corresponds to the intermediate relative momenta p'. The 
deuteron binding energy is represented by Ed, and m stand for the nucleon mass. 

Using the relations for Jacobi momenta given in Eqs. (2.18) and (2.20) the evaluation of 
the permutations in Eq. (2.25) leads to 



{p'm' s2 m' s3 T^\ (qm sl T! \P\ q m° sl T?) \p" m" s2 m" sZ T 2 T'. 



ji a ii ii "\ 



i (p'l (q |qo) |p") 2 ( m ^ m/ s2 m/ S 3 r l r 2 r 3 



+i (p'l (q |qo) |p") 3 {^lm^m'^nr^ 



// I/O ii ii 



+ i ;p 



2 qo 



1 „ 3 



2 qo 



2P + ^qoy i , \ ///l,///l :i r:. 2 , ; '. 

T^P" ~ I qo / 1 1 ( TO sl m s2"4 T l T 2>3 



// 11 11 11 



5 (p' -ir)6 (p" - 7r') ^im^^m^mO^m^m^/TiT^'^rO^^' 

+5 (p + tt) 5 (p" + tt') 5 msim » 2 <5 m ' s2m » 3 (5 m / 3m o i 5 n7 .//^ r »(5 r / r o, 



(2.27) 



where 



tt = -q + qo 



7T = -q- -q . 



(2.28) 



As an aside, the variables are arranged such that each delta function only contains one 
integration variable. Inserting Eq. (2.27) into U^^p, q) in Eq. (2.25) leads to 

^ 1} (P,q)= E (pm s2 m S 3r 2 r3|T(^ p )|7rm°X3^>(^X3^i^i 

m 's3 T 3 

+ E (p^2m s3 r 2 T 3 \T(E P )\ - 7rm' s2 m sl r^) (-n'm^m'^n^ | ^f d ) 



m 's2 T 2 
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Y, (pm s2 m s3 T 2 r 3 \T(E p )\nm sl m' s ry) (n' m'^r' n \ ^ d 

m> s T> 

+ (pm s2 m s3 r 2 r 3 \T{E P )P 23 \ Tzm^m'^r') {it 'm' s m sl t 'n \P 23 l \ ^f d )| 
2 (pm S 2m s3 r 2 r3 |T( J B P )(1 - P 23 )\ nm° sl m' s T°T') l-K'm s m sl T'Ti | 



= £ (pm a2 m a3 r 2 r3 ^m^m'^r') (n'm'^r'n \ < d ) . (2.29) 

i i a ' 
m' s T' 

In arriving at the last equality of Eq. (2.29) we made use of the antisymmetry of the deuteron 
state, ^f d \ and defined a (pm s2 m s3 T 2 T 3 \T(E p ) \ 7rm° 1 m / s r 1 °r / ) a as 



pm s2 m s3 T 2 T 3 \T(E P )\ Tzm^m'j^'j ^ = [pm s2 m s3 T 2 T 3 \T(E p )(l - P 23 )\ nm^m'^r'j . 

(2.30) 

We denote the matrix element a (pm s2 m s3 r 2 T 3 \T(E p ) \ TTm^ 1 m' s T^T') a as physical representa- 
tion of the NN T-matrix, physical meaning that the NN basis states \pm s2 m s3 r 2 T 3 ) contain 
the individual spins and isospins of the nucleons. 

Since the deuteron contains only two definite angular momentum states, it is reasonable 
to apply the standard partial wave expansion 



,M d 

d 



^d d ) = Yl J d PP 2 P'(ls)jm;t)(p(ls)jm;t\^ 1 c 

Isjmt 

= E / d P'p' 2 \p'(U)lM d - 0) MP')- (2.31) 



Here \p'(ls)jm;t) is the standard partial wave basis, and ipi(p') represent the standard s— 
and rf-waves of the deuteron. The projection (7r'm' s m s ir'ri \ ^f^ d ^ on the deuteron state in 
Eq. (2.29) is then explicitly worked out as 

n'm'^r'n | *f d ) = £ / dp' p' 2 (Tr'm'^r'n \p'(ll)lM d ; 0) frtf) 

i J 

= (r / r 1 \0}J2C(lll;^M d -fi) 
in 

x J dp'p' 2 (tt' \p'l/i) (m' s m sl \l,M d -n) ^,(p') 
= c(^0;rV 1 )EC(/ll; / i,M d -/i) 

xY ¥ (tt')C Q^l;m' s m s i^ <5<+m 3 i,M d - / /0*( 7r ') 

= c G^ rV ') c (M i;m >") 

x E C (/ll; M d - m' s - m sl ,m' s + m sl ) 
i 

xY^^^Mn'). (2.32) 
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Finally, inserting Eq. (2.32) into Eq. (2.29) we obtain the first part of the break-up amplitude 
^o (1) (P,q) as 

U^(p, q) = £ a (pm s2 m s3 r 2 r 3 \T(E p )\irm° al m , 9 T°T , ) a C Q|o; r'r^j C m' s m sl 

tti s t 

X C M d - m' s - 771,1, «\ + yj,M d -m' 3 -m s i(7r')V'K 7r ') 

f^ 1 ;^^) (pm s2 m S 3T2r 3 |r(£;p)|7rm° 1 m / s r 1 V 



V2 , V22~ 

X E C ('11; M d " < - "1.1, ™'s + ™-l) ^,M d -<-m sl (7r'M(A (2-33) 
Z 

Eq. (2.33) may serve as a starting point for further expressions for (p, q) to be used in 
the explicit calculations. It shows how Uq 1 \p, q) depends on the PW projected components 
of the deuteron and on the NN T-matrix in a physical representation. In our calculation 
of the Nd break-up process we employ the NN t-matrix in the momentum-helicity basis 
|p; pSA; t) 7Ta [7], where S, t, A are the total spin of the 2N system, the total isospin and the 
helicity. The label na means that the basis state has a definite parity rj^ and is antisym- 
metrized. The connection of the T-matrix elements {pm s2 m s3 r 2 T 3 \T(E P )\ 7rm° 1 m' s r{ ) r / ) a 
to those in the momentum-helicity basis, namely T A (p, 7r; E p ), is given in Ref. [7]. Here 
we want to be more general by letting the nucleon types r 2 , r 3 , r°, T\ being arbitrary but 
employing Kronecker symbols to ensure charge conservation, 

a(T 2 T 3 m s2 m s3 p\T(E p )\T°, -nm^m',*) a 

= JWf-n^- AW E (1 - V,(-) S+t ) C (lit; r 2 r 3 ) C r°, - n ) 

xC Q-S;m s2 m s3 A ) C (- -5; m>x) 

xEW^W-fe'r;^ (2-34) 

AA' 

In the above expression d A , A (9) is a rotation matrix [13]. Using Eq. (2.34) we obtain 
Uq 1 \p, q) in terms of the NN T-matrix elements T AA f(p,7r;E p ) in the momentum-helicity 
basis as 

t/o (1) (p,q) = (M 1;mXl 

x^C(Ill;M d -m^- m sl ,m' s + m sl ) yi,Af d - m ' < - msl (7r')^i( 7r ') 
z 

xg (!-,,(-)-) C (Ii t;T2T3 ) C (ii t ;xf,- Tl ) 
xC (---5;m s2 m s3 A ) C (—5; m>|A;) 

x E <a(^)4a'(^)^( P , tt; £ p ) . (2.35) 

AA' 
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Now let us concentrate on the NN T-matrix elements T]ffi(p, p'; E p ) in the momentum- 
helicity basis. For the calculation of NN scattering it is convenient to choose the z-axis as 
the direction of the initial momenta p'. It is shown in Ref. [7] that in this case the azimuthal 
dependencies of the T-matrix elements can be separated as 

TJ^t (p,p'z; E p ) = e M *Tffi(p,p', cos9; E p ), (2.36) 

which then allowed to reduce the LSE's for Tjf,* (p,p'z; E p ) to be the ones for 
T]ifj{p, p', cos 9; E p ). Thus, in order to calculate Uq 1 ^ (p, q) we have to find a relation between 
7aa*(P) Pj Ep) w hh arbitrary p' and TJ[^f(p,p', cos 6"; E p ), where 9" now depends on p and 
p'. This is done in the following way. First we rotate TJ^f(p,p';E p ) so that p' points in 
the z-direction and then apply the relation given in Eq. (2.36). With R(p') being a rotation 
operator working in momentum and spin space, it follows according to Ref. [7] that 

T^(p,p';E p ) = ™ (p;pSA;t\T(E p )\p';p'SA';ty a 

= ™ (p; pSA; t\ T(E p )R(p') \p% zSA'; t)™ 
= ™ (p; pSA; t\ R(p')T(E p ) \p'z; zSA'; t)™ 

= TO (i2t(p')p;p5A;*|T(S J ,)|p'z;z5A';t) ,re . (2.37) 

Here 

|p; pSA; t)™ = \(l- V,(-) S+t ) (1 + V«) \t) |p; pSA) , (2.38) 

and thus, 

tft(p') | P ; pSA; t)™ = \(l- vA-) S+t ) {1 + Vn) \t)R f (p') |p; P^A) . (2.39) 

The action of R^(p') on the state \p;pSA) leads to two successive rotations as 

R\p') |p; pSA) = R(0, -9', -<f)')R((f)90) \pz; zSA) , (2.40) 

and the result is (see Appendix A for the derivation) given by 

i? f (p') |p; pSA) = e iAn R((/)"9"0) \pz; zSA) 

= e iAn \p";p"SA) , (2.41) 

with 

cos 6" = cos 9 cos 9' + sin 9 sin 9' cos(0 - <p') (2.42) 
sin #V " = - cos 9 sin 9' + sin 9 cos 9' cos(0 -<j>')+i sin 9 sin(0 - 0') (2.43) 

where _Df/ A (0#O) are the Wigner D-function [13]. Inserting Eq. (2.41) and then Eq. (2.36) 
into Eq. (2.37) yields 
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TfJ(p, p' ; E p ) = e~ iAn ™ (p"; p"SA; t| T(£ p ) \p'z; zSA'; t) m 
= e-^TfJ(p",p'z;E p ) 

= e'W'-^TZlf&p', cos9"; E p ), (2.45) 
where the exponential factor e l ( A '<P"- An ) i s calculated as 

UA'^'-AQ) = iA'0" E^=--5-PAfA(0^O)-PAfA'(0 / ^ / O) 

E s N= . s e^*-^d 8 NA (e)d 8 NAI (ff) 



d s A , A {9") 



(2.46) 



Returning to Eq. (2.35), by means of the relation given in Eq. (2.45) we arrive at our final 
expression for U^fa, q): 

xJ2C(lll;M d -m' s - m sl ,m' s + m sl ) Y liMd - m > i - mBl (ir')il) l (ir') 
x £ (l " vA~) S+t ) C (Ht; r 2 r 3 ) C Q±f ; 7?, 



xC Q^S; m s2 m s3 A ) C (—5; rr&mX) 
£ 4 ) A(^)<' ) A'(^)e l(A ' '- An) T A t, i (p, tt, cos£'; E p ), (2.47) 



x 

AA 



with 



cos 9' = cos # p cos 9 W + sin # p sin 9 n cos(0 p — T ) (2.48) 
i(A'0'-An) _ Z)jy=-s e tjV ^ p ^dfj A (9 p )dfj A ,(9 7T ) 

d s A , A (9>) 



(2.49) 



III. RELATIVISTIC KINEMATICS IN THE ND BREAK-UP AMPLITUDE 

In the previous section the break-up operator {7o(p, q) is derived within the framework of 
the nonrelativistic Faddeev theory. Since our goal is to study break-up reactions at interme- 
diate energies, we want to consider the influence of relativistic kinematics. This means, that 
we not only have to employ relativistic energy-momentum relations, but more importantly 
have to reevaluate the Jacobi momenta, carry out corresponding Lorentz transformations 
to the two- and three-particle cm. subsystems, and employ a relativistic description of the 
cross section. For our derivation we adopt the formulation given in Ref. [14]. For clarity we 
will describe the most important steps in detail. 
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A. Jacobi Momenta 



Let a system be described by the energy and momentum vector (E, k) in one frame. 
Then, in a different frame moving with relative velocity u, the system is described by 
(£", k'), connected by a Lorentz transformation L(u) to the first frame, 

(E',k') = L(u)(E,k) (3.1) 

k' = k+ (7- l)(k- u)u-7^u (3.2) 

E' = 7 (£ - k • u) (3.3) 

7 = ~j^= 2 - (3-4) 



vT 



u 



Using these relations we can bring our 3N system from the laboratory frame to the cm. frame 
and find the corresponding Jacobi momenta p and q in the final state and q in the initial 
state. 

As in the previous section we choose without loss of generality the two-nucleon subsystem 
to consist of nucleons 2 and 3, and let nucleon 1 be the spectator. To derive the relative 
momentum p of the subsystem in its cm. frame, we employ the Lorentz transformation 

u = £ (3.5) 

(E' 2 ,p) = (E' 2 X) = L(u)(E 2 M) (3.6) 
(E' 2 , -p) ee (E' 3 ,k' 3 ) = L(u)(E 3 ,k 3 )- (3.7) 

Let us define the following quantities for the 23-subsystem, 

k 23 EE k 2 + k 3 (3.8) 

E 23 ee E 2 + E$, (3-9) 

which are connected by a Lorentz invariant relation as 

E\ 3 - k 2 23 ee M 2 3 > 4m 2 . (3.10) 

Here M 23 is called the invariant mass of the 23-subsystem and equals the total energy of 
the 23-subsystem in its cm. frame. According to Eq. (3.2) and the transformations given 
in Eqs. (3.5)-(3.7), the Jacobi momentum p is given as [15] 

p .I( k2 - k3 )-I k23 (|^|L). (3.11) 

The last term in Eq. (3.11) exhibits the relativistic effect in the definition of the relative 
momentum. Since in the cm. frame the energy E 2 is equal to E' 3 , the total energy M 23 in 
the 23-frame is given as 



M 23 = E' 2 + E' 3 = 2E' 2 = + k' 2 2 = 2^m 2 + p 2 . (3.12) 

Starting with Eq. (3.12), employing energy and momentum conservation together with 
Eqs. (3.8)-(3.10) the magnitude of p is calculated as 
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P=\\l™ 2 d- 2m2 + ^a{E lab - E x ) - 2E lab E 1 + 1k lah k x cos6 lab . (3.13) 

Knowing the total energy M 23 in the 23-frame one can calculate the kinetic energy Ep in 
the 23-subsystem, that is 

E p = M 23 - 2m = 2^/m 2 + p 2 - 2m. (3.14) 

Thus the NN T-matrix elements T][[t(p, n, cos 9'; E p ) in Eq. (2.47) will be calculated for the 
energy E p given in Eq. (3.14). 

In order to obtain the Jacobi momentum q we apply the following Lorentz transformation 

u=^ (3.15) 

Eq 

(E' 1 ,q) = (^,k , 1 )=L(u)(S 1 ,k 1 ) (3.16) 
(E' 23 , -q) = (E' 23 , k' 23 ) = L(u)(E 23 , k 23 ), (3.17) 

where E is the total energy in the laboratory frame, 

E = m d + E lab = E 1 + E 2S . (3.18) 

Similar to p, the Jacobi momentum q acquires a relativistic correction term and is given by 

q=-(k 1 -k 23 ) + ^^ Eq + Mq (^-^ 23 )j. (3.19) 

Here M is cm. total energy or the invariant mass of the 3N system 

M = E' 1 + E' 23 , (3.20) 

which is connected to the laboratory total energy E and the laboratory total momentum 
k; a 6 by the following Lorentz invariant relation 

El - kl b = Ml > 9m 2 . (3.21) 

The energies E[ and E' 23 are given in terms of the magnitudes of Jacobi momenta p and q 

as 



E[ = ^m 2 + k' 1 2 = ^m 2 + q 2 (3.22) 

E' 23 = ^M 2 3 + k' 23 2 = ^4(m 2 +p 2 )+q 2 . (3.23) 

With Eqs. (3.16)-(3.18), (3.20) and (3.22)-(3.23) we get for the scalar product in Eq. (3.19) 

(ki - k 23 ) • k tab = k\- k\ 3 

= El - E\ 3 - m 2 + M| 3 

= (Ei — E 2 3)E — (E[ — E 23 ) 

= (Ei - E 23 )E - (E[ - E' 23 )M . (3.24) 



12 



Inserting this result into Eq. (3.19) leads to a final expression for q, 

H 2 V i! 2M V E + M V i) ) 

- -(kx - k 23 ) - -k, ab ^ e^M, ) 

For the last equality total momentum conservation was employed. With the help of 
Eqs. (3.22) and (3.23) we can write M as 

M = E[ + E' 23 = ym 2 + q 2 + yj 4(m 2 + p 2 ) + q 2 . (3.26) 
Starting from Eq. (3.26) we obtain after some algebra the magnitude of q as 

'{M^ - (5m 2 + V)} 2 - 16m 2 (m 2 + p 2 ), (3.27) 



* 2M 

where the value of M can be calculated from Eq. (3.21) as 

Mo = \[E 2 - kf ab = yjm 2 + m 2 + 2m d E lab . (3.28) 

Finally we turn to the initial state, calculate the Jacobi momentum q , the energies E[ ab 
of the incoming nucleon and E' d of the deuteron in the cm. frame. Replacing in Eq. (3.19) the 
quantities ki with k; ab , E 1 with E tab , k 2 3 with zero and E 2 3 with m^, the Jacobi momentum 
qo is given as 

q ° = l klab + 2lt {eV+M " {Elab ~ md) ) 
= ±k lab + ^±(E -M - (E lab - m d )) 

= ^ la , (3.29) 
The energies E' lab and E' d are obtained using Eq. (3.3) as 

£ " = t( £ "-t) = ^ (m2+mAt) (330) 

E 1 

E ' d = W md = M~J KElab + m ^ md - ^ 3 - 31 ^ 

It can be shown that E[ ab and E' d sum up to M as required by total energy conservation in 
the cm. frame. 
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B. S-Matrix 



For a correct treatment of relativistic kinematics, we still need to connect the S-matrix 
element in the laboratory frame to the Nd break-up amplitude. Suppressing all discrete 
quantum numbers we define the S-matrix element for the break-up process in the laboratory 
frame as 



S , (ki,k 2 ,k 3 ) = (kik 2 k 3 \S\ki ah k d ) 



and in the cm. frame as 



S(p,q) = (pq|S|qo) = (pq|S|qo>- 



(3.32) 



(3.33) 



In the above equation we omitted the deuteron state \^d) m the initial state, but showed 
its momentum k^ in Eq. (3.32), though its value is zero. The momenta p, q, and qo are 
the Jacobi momenta calculated in Eqs. (3.11), (3.25), and (3.29). The states |kxk 2 k 3 ) and 
\ki ab k d ) are related to the states |pq) and |q ) by the following relations, 



Ikxkaka) = IkOJkaks) 

= J-^ka^lkx) |pk 23 ) 
= J _ 5(k 2 ,k 3 )J"^(ki,k 2 3)|pq)]ki + k 23 ) 
\k lab k d ) = J~2(k Za6 ,k d )|qo)|k Za6 ). 



(3.34) 
(3.35) 



The Jacobian J(k 2 , k 3 ) of the transformation from the variables (k 2 , k 3 ) to (p, k 23 ) is given 
by [14] 



</(k 2 ,k 3 ) 



<9(k 2 ,k 3 



d(p, k 



2: j ,j 



E 2 E 3 M; 



2.3 



E 23 E' 2 E 3 



AE 2 E 3 
E 23 M 23 



(3.36) 



where the last equality results by means of Eq. (3.12) for M 23 . Similarly, the Jacobians 
J(k 1; k 23 ) in Eq. (3.34) and J(ki ab ,k d ) in Eq. (3.35) are given as 



J(ki,k 23 ) 

J (klab, k d ) 



d(k u k 23 ) 



<9(q, kx + k 23 ) 
d(k lab , k d ) 



9(q , k iab ) 



E\E 23 M 
Eq E[E 23 
Ei ab m d M 
E o E' lab E' d 



(3.37) 
(3.38) 



Thus, one can finally relate S(ki, k 2 , k 3 ) to ^(p, q) as 

<S , (ki,k 2 ,k 3 ) = (kik 2 k 3 |S , |ki a6 k d ) 

= J~5(k 2 ,k 3 )J-5(k 1 ,k 23 )J-5(k Za6 ,k ( i)(k 1 + k 23 |(pq| 1 S|qo)|k te6 ) 
= 5(ki + k 23 - ki ab ) { J(k 2 , k 3 ) J(ki, k 23 )J(k lab , k d )y^ S(p, q), 



(3.39) 



where the delta function ensures total momentum conservation. 

As last step we need to connect S(ki, k 2 , k 3 ) to the break-up amplitude Uo(p, q) defined 
in Eq. (2.12), and which is related to S(p, q) as 
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S(p, q) = -2m5(E[ + E' 23 - E[ ab - E' d )U (p, q). (3.40) 

Inserting Eq. (3.40) into Eq. (3.39) this gives 

5(k 1; k 2 , kg) = -2m5(k 1 + k 23 - hab)S(E[ + E' 23 - E' lah - E' d ) 

x { J(k 2 , k 3 ) J(k!, k 23 ) J(k ia5 , k d )}^ C/ (p, q). (3.41) 

By means of the identity 

5(k x + k 23 - k^)*^ + E' 23 - E' lab - E' d ) 

= 5(ki + k 23 - kfafc)^ + £ 23 - ^ - m d ) f 1 ± f 23 ± ^ ± ^ , (3.42) 

iii + -fc 23 + biab + m d 

we obtain the relation between S(ki, k 2 , k 3 ) and U (p, q) as 

S(ki, k 2 , k 3 ) = -27ri5(k! + k 23 - k ia6 )5(£i + E 23 - E lab - m d ) 

x E' l + E' 23 + Ej ab + E' d i 

= -2 7 r^(k 1 + k 23 - k lab )5(E 1 + £ 23 - E lab - m d )r(p, q)U (p, q). (3.43) 
Here the function T(p, q) is defined as 

r( P ,q) = ftf't E } h t E ' d {J(k 2 ,k 3 )J(k 1; k 23 )J(k, afci k d )}^ 
£/! + £/ 23 + £/ iafe + m d 



\ M 2 ?,E' l E' 23 E' lab E' d 
" V ±E x E 2 E 3 E lab m d ' 

In the nonrelativistic limit the function T(p, q) equals 1. 



(3.44) 



IV. OBSERVABLES IN THE PROTON-NEUTRON CHARGE EXCHANGE 

REACTION 

So far we derived the break-up amplitude U (p, q) in first order in its most general form. 
As application of a proton-deuteron break-up process we consider the (p,n) charge exchange 
reaction. In the experiments we are going to analyze only the scattered neutron is detected. 
Thus, when calculating the observables of this reaction, all possible directions of the two 
undetected protons must be taken into account. This is accomplished numerically by inte- 
grating over the relative direction p between the two protons. In our numerical application, 
we consider the spin averaged differential cross section in the d(p, n)pp reaction and selected 
spin observables. These are the neutron polarization Pq in the d(p, n)pp reaction, the an- 
alyzing power A y in the d(p, n)pp reaction and the polarization-transfer coefficients Dij in 
the d(p, n)pp reaction. Comprehensive descriptions and derivations of these observables can 
be found in e.g. Ref. [2]. 

The nonrelativistic cross section is given as 



d 5 a _ ^ Hfpki 1 
(■//•.; <7k i -h'1,,1, 6 „, .. 



= (2-) 4 ^Z E W(p,q)| 2 , (4.1) 
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where ki determines q via Eq. (2.8) and p via Eqs. (2.26) and (2.8). The relativistic cross 
section is given by 

da = ( 27r ) 4 ^5(k 1 + k 23 - \^i ab )5{E l + E 23 - E lah - m d ) 

Rab 

xr 2 (p,q)|f/ (p,q)| 2 rfk 1 rfk 2 rfk 3 , (4.2) 

where the S-matrix element from Eq. (3.43) is used, and all energies obey the relativistic 
energy-momentum relation. The cross section for the inclusive Nd break-up process is cal- 
culated as function of the direction kx and the kinetic energy E^i — Ei — m of the detected 
nucleon. Using the relation 

dki = dhk^dki = E^dE^ = E^dE^d^ (4.3) 

leads to 

dcr , n A Ei ah Eiki 



(2tt) z 

dExdkx k lab 

x J dk 2 dkg<S(ki + k 23 - k lab )5(E 1 + E 23 - E lab - m d )T 2 (p, q)\U (p, q 

= (2tt) 



uEi ab Eiki 



hab 

x J dk 23 dpj(k 2 , k 3 )<5(k 1 + k 23 - k ia5 )5(E 1 + E 23 - E lab - m d )r 2 (p, q)|£/ (p, q) 
^yE^M I d^^^s^ + E^-E^-m^^^lUoip,^ 

Klah J &1?,Mo 3 



= W I m I rfp^3r 2 (p,q)|f/o(p,q)| 2 . (4.4) 
In addition we used Eq. (3.36) for J(k 2 ,k 3 ) together with 

dE 23 = d^Ml 3 + k\ 3 = d^A(m 2 +p 2 ) + k 2 23 = ^-dp. (4.5) 

-^23 

Defining a function p(p, q) by 

s 2E^^ = fiAS"^ , (4.6) 

M 23 zm^ 

allows to write the relativistic cross section similar to the nonrelativistic one, that is 

^ -M 4 %^i £ /#l«P,q)P. (4.7) 



rf_Eic?ki 



Here we restore the summation over final spins and the averaging over initial spins states. 
In the nonrelativistic case, the function p(p,q) reduces to m 3 , leading to Eq. (4.1). Next 
we give the polarization P , the analyzing power A y and polarization-transfer coefficients 

£>y = 4rr{Ub(^-j)^-i)}as 
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Po = 


±Tr{u Ul[ 


>•*)} 




(4.8) 


Ay = 


w Tr i Uoitr 


■ml} 




(4.9) 


D n 'n 


^-Tr{U (<r 


■h)Ut(<T 


■n)} 


(4.10) 


Dg's = 


^-Tr{u (a- 




•s')} 


(4.11) 


D Vs = 




■s)Uo(<r 


•1')} 


(4.12) 


Dsn = 




•i)[/ V 


■0} 


(4.13) 


Dm = 




• i)tf}(«7 • 




(4.14) 



where 



l -Tr{u Q Ul). (4.15) 



For simplicity we suppressed the p-integration in Eqs. (4.8)-(4.15). The unit vectors n, 1, s, 
n', 1' and s' are defined as 

n = n l = * labX *\ , i = k, a6 , SEnxl, i' = k 1; s' = n' x 1'. (4.16) 

\klab X Kl| 

For the explicit calculation of the Nd break-up amplitude Uo(p, q), Eq. (2.47), we need 
the NN T-matrix elements. They are obtained by solving the Lippmann- Schwinger (LS) 
equations for a given NN potential in 3D as described in Ref. [7] at fixed momenta, angles, 
and energies. The matrix elements TJ[^f{p, ir, cosO'; E p ) are then obtained by interpolating 
in 71, cos#', and E p . This is more economic than solving the LS equation every time for 
the corresponding energies and initial momenta. Similarly the partial wave deuteron wave 
function components are calculated once and interpolated to the momenta ir'. For all the 
interpolations we use the modified cubic Hermite splines of Ref. [16]. Typical grid sizes for 
our calculations are 40 points for E p (and p), 80 points for cos8', and 50 points for ir. This 
grid is then used for all our break-up calculations from 5 to 500 MeV. 

V. RESULTS AND DISCUSSION 

Our numerical calculations are based two different NN potentials, the Bonn-B [17] and 
AV18 [18] potentials. In the discussion of our results we consider three different aspects. 
First, since we present a new way of calculating the break-up process, namely without partial 
wave decomposition, we need to compare our results to those obtained in traditional partial 
wave based calculations. After having established the feasibility and correctness of our 
calculations, we will address two physical questions, namely the importance of rescattering 
contributions in the (p,n) charge exchange reaction at a moderately high energy and the 
effect of different descriptions of the kinematics, i.e. we compare the relativistic treatment 
introduced in Section III to the nonrelativistic description. 
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A. Comparison with Partial- Wave Calculations 



In this section we compare our 3D calculations with traditional PW calculations at dif- 
ferent energies. The aim here is twofold. First, we want to convince ourselves that our newly- 
developed 3D formulations agrees with well established PW calculations [19]. Secondly, we 
want to find out from which energy regime on our 3D method surpasses the state-of-the-art 
PW calculations. As already mentioned in the introduction, in PW calculations the number 
of partial waves necessary for a converged result proliferated as the energies increases, leading 
to limitations with respect to computational feasibility and accuracy. Our new 3D scheme 
does not suffer from these limitations, the computational effort is in principle independent 
of the projectile energy considered. 

Our first comparison is carried out at a proton incident energy of Ei a b = 100 MeV and 
a neutron laboratory scattering angle Q\ a b = 13°. The calculations, which are based on 
the Bonn-B potential are given in Fig. 1, which shows the differential cross section, the 
analyzing power A y and the polarization-transfer coefficient D s i. The solid lines represent 
our 3D calculations, the dashed lines the corresponding PW calculations [19]. Here we use 
the notation j for the highest 2N total angular momentum taken into account in the PW 
calculation, and J for the highest 3N total angular momentum. The figure shows that 
both lines are almost indistinguishable, thus validating our new scheme. At this point we 
also would like to mention, that the channels used in the PW calculation, namely j = 7 
and J = 31/2 constitute todays limits for a PW calculation. In addition, we carried out 
comparisons at lower energies, e.g. at 16 MeV, where a PW calculation with j = 5 and 
J = 31/2 is in perfect agreement with our 3D calculations. 

Next, we turn to a slightly higher projectile energy, E iab = 197 MeV, and carry out 
the same comparison. The results for the differential cross section, the analyzing power 
A y and the polarization-transfer coefficient D s i are shown in Fig. 2, where the solid line 
represents our 3D calculation. The PW calculations are shown with increasing number of 
partial waves from j = 5, J = 25/2 to j = 7, J = 31/2. The peak of the differential cross 
section reveals that each additional angular momentum of the PW calculation results in an 
additive contribution, but even the highest possible number deviates about 9% from our 
3D result. This is the most extreme case, for the analyzing power A y and the polarization- 
transfer coefficient D s i the PW calculation with j = 7 and J = 31/2 agrees reasonable well 
with our 3D result. It is interesting to note that for D s i the 2N total momentum j is much 
more important to reach convergence than the total 3N momentum J . 

At this point it is appropriate to make some general remarks. In this work we restrict 
our 3D approach to the leading term in the Faddeev multiple scattering series. Thus we 
have no insight whether the fully summed series would lead to a better agreement of 3D 
and PW approach. We also restrict ourselves to semi-exclusive processes, and can not draw 
any conclusions on 3D and PW calculations with respect to elastic scattering observables or 
full break-up observables. But nevertheless, the comparison in Fig. 2 indicates that around 
200 MeV the convergence of the PW approach has to be carefully checked. 
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B. Contributions from the Rescattering Terms 



One of the arguments to study the semi-exclusive Nd break-up process in first order in the 
NN T-matrix is that at higher energies the rescattering term generated by the solution of the 
full Faddeev equations become less important. For a comprehensive study of the importance 
of those rescattering terms it would be necessary to compare first order calculations with 
full Faddeev calculations over a wide range of projectile energies and for many different 
experimental situations. Unfortunately we can not do this at the present stage, since three- 
dimensional full Faddeev calculations do not yet exist, and traditional, partial wave based 
Faddeev calculations are limited in their energy range. Thus we take as compromise a 
medium energy of about 200 MeV, and compare the (p,n) charge exchange observables 
calculated in our first order 3D approach with the ones obtained from a full, partial wave 
based Faddeev calculation. We choose the proton energy Ei a b=197 MeV, since there exist 
recent measurements [8]. 

Our calculations are based on two different NN potential models, namely Bonn-B [17] 
and AV18 [18]. Both potentials are defined below 350 MeV nucleon laboratory energy, which 
corresponds to a NN cm. energy of 175 MeV. In the Nd break-up process in first order the 
NN cm. energy available to the two-nucleon subsystem is fixed in terms of the laboratory 
momentum of the final nucleon and the projectile energy. For a projectile energy of about 
200 MeV in the pd scattering process, the maximum NN cm. energy in the two-body 
subsystem is about 133 MeV. Thus, our calculations employs the NN models in an energy 
regime, where they are perfectly well defined. Of course, the two potential models exhibit 
differences in the description of the NN phase shifts. The model AV18 is one of the socalled 
high-precision potentials, describing the NN data base with a x 2 /datum fa 1, whereas Bonn- 
B has a slightly higher x 2 / datum value. Thus, there are on-shell differences between those 
two models, which should lead to differences in the Nd break-up observables. 

In Figs. 3 and 4 we compare the 3D calculations with PW based full Faddeev calculations 
[19] at 197 MeV proton energy. We show the spin averaged differential cross section, the 
analyzing power and polarization-transfer coefficients at two angles, Q\ a b = 24° and 6^ = 37° 
together with experimental results from Ref. [8]. The PW full Faddeev calculations use 
j = 5, J = 31/2 for the AV18 NN potential, and j = 4, J = 31/2 for the Bonn-B NN 
potential. Since the solution of the PW full Faddeev equations is more involved we restrict 
ourselves to a lower number of partial waves. However, the number of partial waves is 
sufficient to study the importance of rescattering terms in the multiple scattering series at 
this energy. The first obvious difference between the two calculations is the appearance of the 
final state peak in the differential cross section, which of course is solely due to rescattering. 
Furthermore, we see that rescattering contributions have the general tendency to push the 
peak of the differential cross section down, though the size of the effect depends on the 
angle and the potential. However, the rescattering terms do not affect the position of the 
peak. We see that the peak is shifted further away from the data the larger the neutron 
scattering angle becomes. For both angles, the analyzing power A y shows the largest effect 
of rescattering for small neutron energies, which can be expected, since interactions between 
outgoing particles should be larger, when their relative energy is smaller. In both cases, A y 
can only be satisfactorily described when rescattering terms are taken into account. For the 
polarization-transfer coefficient D u at 9 iab = 24° the situation is similar, rescattering effects 
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are largest for small neutron energies. For D ss at Q\ a b = 37° none of the calculations is able 
to capture the general shape of the data, rescattering effects are visible, but they do not 
affect the general shape of the curve as this is the case in the other observables shown in 
Figs. 3 and 4. From these consideration we have to conclude, that at a projectile energy 
Ei ab ~ 200 MeV rescattering terms still give considerable contributions to the full pd break- 
up amplitude and hence cannot be neglected. Due to the lack of calculations based on the 
full Faddeev equations at higher energies, we can not carry out corresponding studies at 
higher energies. 

C. Relativistic Effects 

In this section we study the effects of relativistic kinematics in the break-up amplitude 
and follow the formulation derived in Section III. We also want to take full advantage of 
our 3D formulation and carry out calculations at proton incident energies higher than 200 
MeV, a regime where partial wave based Faddeev calculations become less competitive. Of 
course we also realize that the NN potentials from which our NN T-matrix is obtained are 
strictly speaking out of their range of validity, i.e. they do not include important Delta 
degrees of freedom. A NN laboratory energy of 350 MeV roughly corresponds to a proton 
incident laboratory energy of 260 MeV in the pd scattering process. A comparison of the NN 
scattering observables with the calculated ones shows that even at NN laboratory energies 
higher than 350 MeV the agreement with data deteriorates relatively slowly. Nevertheless, 
this can lead to deficiencies in describing the Nd break-up process at Ei ab >260 MeV. 

At first we investigate the effect of relativistic kinematics on the break-up observables 
at a low energy, = 100 MeV, where it is expected to be small. In Fig. 5 we show the 
differential cross section and A y at E\ a b = 100 MeV and a neutron laboratory scattering 
angle Qi a b = 24°. In the cross section effects are only visible in the quasi- free-scattering 
(QFS) peak, but in general one can say that around 100 MeV relativistic effects are small, 
and certainly not the dominant correction to worry about. 

Going to a higher energy, = 197 MeV, the relativistic effects increase considerably. 
In Figs. 6 and 7 we show the cross section and A y and two polarization-transfer coefficients 
at neutron laboratory scattering angles d\ a b = 24° and d\ a b = 37°. Here the QFS peak is 
visibly enhanced by the use of relativistic kinematics. More importantly, it location is shifted 
towards smaller neutron energies, and is now in better agreement with the experimentally 
determined peak location. As far as the spin observables are concerned, the relativistic 
corrections show the largest effect for the higher neutron energies. 

A proton energy of E\ a b = 346 MeV is the next higher energy at which the (p,n) charge 
exchange reaction is measured [9]. In Fig. 8 we display 3D calculations with nonrelativistic 
and relativistic kinematics for a neutron scattering angle 8i ab = 22°. Again we observe an 
increase in the QFS peak and a shift to lower neutron energies. Here we would like to 
point out that there is an uncertainty in the data as far as the location of the QFS peak is 
concerned. In the experiment there is an uncertainty of the energy, at which the pd break-up 
process exactly occurs. For example, due to the thickness of the target the proton may have 
lost some of its energy before it hits and breaks the deuteron apart [19]. In this case, the 
break-up the process occurs at an energy slightly different from the calculated one. The 
effect of relativistic kinematics on the spin observables is clearly more pronounced compared 
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to E tab = 197 MeV. We also calculate the break-up process at E iab = 495 MeV, though here 
the uncertainty with respect to our input NN interactions is largest. In Fig. 9 we show the 
differential cross section, A y and D a at E iab = 495 MeV for a neutron laboratory scattering 
angle 9i ab = 18° together with the measurements from Ref. [21]. Here we clearly see that the 
corrections due to relativistic kinematics push the QFS peak towards lower neutron energies, 
and the location of our calculated peak agrees with the measured one. The effects on the 
spin observables are now also quite sizable. 

With this study we can qualitatively indicate that relativistic effects become important 
when going to higher energies. However, we can not make any definite statements, since we 
only consider relativistic kinematics. We have not considered effects resulting from boosting 
the NN T-matrix [10] and Wigner rotations of the spin [11]. Those effects in principle could 
counterbalance the kinematic effects. It also remains to be seen, how important rescattering 
effects will be at those higher energies. Our calculations based on the first order term and 
relativistic kinematics overestimate the differential cross section. That could imply, that 
rescattering still plays an important role at those energies. 

VI. SUMMARY AND CONCLUSIONS 

We formulate and calculate the Nd break-up process based on the Faddeev scheme in 
first order in the multiple scattering expansion in a three-dimensional fashion which does not 
rely on any partial wave decomposition. The leading term for the Nd break-up amplitude is 
derived in a representation, which uses directly the momentum vectors. This representation 
can be connected to the momentum-helicity basis, in which we solve for the NN T-matrix in 
a 3D fashion. Special care has to be taken when rotating the NN T-matrix elements, which 
occur with arbitrarily oriented initial momenta in the Nd break-up amplitude, such that the 
NN initial relative momenta point into a fixed z-direction. This is needed since two nucleon 
LS equation for the NN T-matrix is solved in a basis where the arbitrary z-axis points into 
the direction of the initial momenta. This leads to an intricate additional phase factor. 

As specific application of our new formulation we calculate the (p,n) charge exchange re- 
action in the proton-deuteron break-up process. Here only the outgoing neutron is detected 
after the break-up. Our calculations concentrate on spin averaged differential cross sec- 
tions, neutron polarizations, proton analyzing powers, and polarization-transfer coefficients 
at different energies. 

First we carry out calculations of observables for the leading order term in the NN 
T-matrix at energies which are accessible to traditional, partial wave based Faddeev calcu- 
lations. The aim here is twofold. First, we need to establish the numerical accuracy and 
feasibility of our new formulation. We establish both by comparing observables calculated 
in both schemes at a proton incident energy E iab — 100 MeV, where we find excellent agree- 
ment between both calculations. At E[ ab = 200 MeV we find some slight deviations between 
the two schemes, especially in the quasi-free peak of the differential cross section. This can 
be identified as the onset of a lack of convergence using the typical and feasible number of 
partial waves in the traditional partial wave based calculation in that particular observable. 
From that we conclude that starting at about 200 MeV, the convergence of partial wave 
based calculations has to be checked carefully. 
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Second, we want to investigate the importance of rescattering terms at a moderately 
high energy. Of course, we need full Faddeev calculations here. Since those do not yet exist 
in a 3D formulation using realistic NN potentials, we have to resort to partial wave based 
full Faddeev calculations. This of course limits the energy regime we can study. Thus, 
we compare our calculations at Ei ab = 197 MeV to the PW full Faddeev calculations. We 
find that at this energy rescattering effects are still important, and are mostly visible in 
the cross section and the analyzing power. In addition, we find that the PW full Faddeev 
calculations provide a reasonable description of the (p,n) charge exchange reaction at 200 
MeV. However, we also can detect one obvious deficiency in both schemes, at larger neutron 
laboratory scattering angles the QFS peak is located at slightly too high neutron energies 
compared to the data. 

This leads to the next topic we investigate, namely the effect of relativistic kinematics in 
the Nd break-up reaction. Here we have to employ not only relativistic energy-momentum 
relations, but also need to reevaluate the Jacobi momenta by carrying out corresponding 
Lorentz transformations to the two- and three-particle cm. subsystems, and employ a 
relativistic description of the cross section. We compare our 3D calculation based on non- 
relativistic kinematics with the corresponding one based on relativistic kinematics. Though 
there are no sizable effects at 100 MeV proton incident energy, we find, that at 200 MeV 
visible effect occur, mainly the differential cross section. Its magnitude increases, but most 
importantly, the QFS peak is shifted to the experimentally determined one using relativistic 
kinematics. Since our calculations are as easily carried out at 300 or 500 MeV as at 200 
MeV, we perform calculations at £^(,=346 MeV and 497 MeV, where experimental data are 
available. We find that at those higher energies the effects due to the relativistic kinematics 
are considerably larger than at 200 MeV. They are now visible not only in the cross sec- 
tion but also in the spin observables. Even at Ei ab =500 MeV this specific feature prevails, 
namely that the QFS peak is shifted to lower neutron energies and coincides now with the 
experimentally determined one. However, its magnitude is larger. With these finding we can 
qualitatively indicate that relativistic effects become increasingly important when consider- 
ing the Nd break-up reaction at higher energies. However, we have to exercise some caution 
in the interpretation of our findings, since we have not considered dynamical relativistic 
effects, such as boosting of the NN T-matrix. It also remains to be seen how important 
rescattering effects will be in the higher energy regime. 

Summarizing, for the first order term in the multiple scattering series in the Faddeev 
scheme the 3D approach has proven to be a viable alternative to the established partial wave 
based calculations. When entering the intermediate energy regime it may be the approach 
having the most promise of being successful in the near future, due to the intrinsic limitations 
with respect to computational feasibility and accuracy faced by partial wave calculations at 
higher energies. It is also clear that the 3D approach, though having a well defined roadmap 
ahead, is still facing extensive development needs. The full Faddeev equations will have to 
be solved, with the inclusion of three-nucleon forces, which may play a more dominant role 
at higher energies. Furthermore, though we consider the effects of relativistic kinematics, we 
have not taken into account the corresponding dynamical effects. And last, but not least, 
the underlying input of any 3N calculation, namely the two-nucleon force, is by far less 
developed at higher energies than it is for energies below the pion production threshold. 
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APPENDIX A: TWO SUCCESSIVE ROTATIONS 



In this appendix we evaluate the rotation of the state \p;pSA) as 

R j {p) |p; PSA) = R{0, -6', -<j>')R(<j>BQ) \pz; zSA) . (Al) 

First, we give a few basic definitions and relations required to follow the calculation. More 
details about rotation can be found in e.g. Ref. [13]. 
A general rotation operator R(p) is defined as 

R(P) = R(<f)60) = e -iJ*<f> e -iJy^ (A2) 

where J z , J y are the z- and y-components of the angular momentum operator J and (6,<f>) 
the rotation angles that determine the direction of p. This operator rotates the angular 
momentum state \zjm) into the state \pjm), 

\pjm) = R{p)\zjm) = J2Di m (p)\zjm'), (A3) 

m' 

where D J m , m (p) are the Wigner D-function defined as 

DLm(p) = D j m , m {m = (Zjm'\R(p)\zjm). (A4) 

A rotation R(a(3 / y) corresponds to a change of the Cartesian coordinates r describing the 
state. The new Cartesian coordinates r' are related to the old ones r as 

r = M(a/37)r, (A5) 

where the rotation matrix M{a(3~f) is given as 

(cos a cos j3 cos 7 — sin a sin 7 sin a cos (5 cos 7 + cos a sin 7 — sin (3 cos 7 \ 
— cos a cos (3 sin 7 — sin a cos 7 — sin a cos (3 sin 7 + cos a cos 7 sin (3 sin 7 
cos a sin (3 sin a sin (3 cos /3 / 

(A6) 

1. Two Successive Rotations in Momentum Space 

We denote the rotation operator in momentum space as Rl(p), which is given in term 
of the orbital angular momentum operator L as 

R L (p) = R L {<f)60) = e -iL,4> e -iL y e_ ( A? ) 

A momentum state |p) with p pointing in the direction (9, (f>) can be expanded in partial 
waves as 

\p} = J2\plm}Y l * m (e,<j ) ), (A8) 

lm 
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where \plm) is defined to be quantized along the z-axis. The state |p) can be obtained by 
rotating a state \pz) as follows, 

R L (P)\pz) = i^0O)£|jrfm)l^(O,O) 

lm 



/2/ + 1 



47T 

2/TT 



= EE \P l ' m ) (%'m\R L (000) | zlO) 

I I'm 



+K 



lm 



121 + 1 

47T 



= EbMW>) 

= |P>- (A9) 
Here we used the relation between the spherical harmonics and the Wigner D-functions, 



i:(m) = ||^w (aio) 

Now we rotate the state |p) with an inverse rotation operator R~l 1 {p') = R L (p') = 
R L (0, -9', -<f)'). It follows that 

4(p» = 4(p')EW^(M) 

lm 

= EE b/W)(z/W|4(p')|z/m)F^(^0) 
= EEl/ m '}(^'m'|4(p)|zM (zZm |p) 

/m I'm' 

= J2\pl'm')(zl'm'\Ri(p')\p) 

I'm' 

= J2\pl'm')(zl'm'\ p") 

I'm' 

= ElA')C(«''f) 

I'm' 

= R L (P")\pz), (All) 
where we have defined a direction p" to be determined by p and p' according to 

| P ") = 4(P')|P>- (M2) 
Inserting Eq. (A9) into Eq. (All) this leads to 

R[(P')R L (P)\pz) = R L (p")\pz). (A13) 

Hence, the two successive rotations R} l {P')Rl{P) applied to the state \pz) can be replaced by 
the single rotation R L (p"). Consequently any number of successive rotations in momentum 
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space can always be replaced by one rotation with the corresponding rotation angles. The 
angles (9", 4>") of p" are determined by the angles (9, 4>) of p and (9', 0') of p' by 

cos 9" = cos 9 cos 9' + sin 9 sin 9' cos(0 - 0') (A14) 
sm9"e ir = -cos9sm9' + sm9cos9' cos((f)- <f)') +ism9sm((f)- (f)'), (A15) 

and are obtained from the rotation matrices of the Cartesian coordinates, which correspond 
to the rotations in Eq. (A13). Such a rotation matrix M(a(3~f) corresponding to R{af3^) is 
given in Eq. (A6). 



2. Two Successive Rotations in Spin Space 

We denote the rotation operator in spin space as Rs(p), which is given in term of the 
total spin operator S as 

R s (p) = R s (<j)90) = e ~ iS ^e- lS ^. (A16) 

The rotation identity given in Eq. (A13) may not apply in spin space. Therefore, we 
evaluate two successive rotations in spin space, independent of the evaluation in momentum 
space. We compare the rotated spin state or the helicity state |p"SA)j with |p"SA)jj given 
by ' ' 

|p"S'A) I = J Rs(p")|25 , A) (A17) 
|p // 5A) n = 4(p / )i?5(p)|z5A). (A18) 

It should be pointed out that here the relation between (9",(f)"), (9',(f)') and (#,0) given in 
Eqs. (A14) and (A15) is still valid, since transformations of the Cartesian coordinates are 
the same in both momentum space and spin space. 

Both states |p"SA)j and |p"SA)jj are eigenstates of the helicity operator S • p" with 
eigenvalue A, as can be shown as follows: 

S • p"|p"SA> T = Rs(i>")S ■ zRl(p")R s (p")\zSA) 
= R s (p")S ■z\zSA) 
= Ai2 s (p")|zS'A) 

= Alp'^A)! (A19) 

S • p"|p"SA> n = 4(PWP)S • z4(p)i? 5 (p')4(p' )Rs(p)\iSA) 
= Rl(p')R s (p)S-z\zSA) 
= ARl(p')R s (p)\zSA) 

= A|p"SA) n . (A20) 

Moreover, because of the unitarity transformations in Eqs. (A17) and (A18) the two states 
have the same norm and can at most differ by a phase factor. 

The helicity states |p"SA)j and |p"SA)jj are expanded in the spin states \zSA) as 



26 



|p"SA)i = R s (p")\zSA) = ^ \zSA')D s A , A ( ( j ) "9"0) 

A' 

\p"SA) u = RUp')Rs(p)\zSA) 

= \zSA')(zSA'\Rl((f)'6'0)\zSN)(zSN\Rs((f)eO)\zSA) 



(A21) 



A'N 



EI^A'>E^A'(^0)Dj A (^0) 



A? 



£|zSA')Xf, A (0V'O), 



(A22) 



where 



Xf, A (0V'O) = E^a'(^'0)^a(^0). 



(A23) 



N 



Therefore, instead of comparing |p // 5 r A)j with Ip'^A) jj we compare Df, A ((J)"6"0) with 
Xl, A ((f)"9"0), since these are known functions. We have two spin cases S = and S = 1. 
For 5 = the spin state is rotationally invariant and thus we can immediately get 

X O ° O (0"£"O) = DM"e"0) = 1 (A24) 

and correspondingly 

|p"00) T = |p"00) n = |z00>. (A25) 
For S — 1 we make use of a symmetry relation for the Wigner D-functions given as 

I%J<*M) = (-) m '- m Dl m ,- m (*f3l), (A26) 

allowing to leave out the case with initial helicity A = — 1 and consider only six cases with 
A' = 1,0,-1 and A = 1,0. 

The Wigner D-function D A , A ((f)0O) is given as 



D\(f>eo) 



For A = it follows that 



p — i(f> 1+cos 9 — icf> sing r —i<t> 1— cos 9 \ 

C 2 V2 2 I 

sing „„„/) sin 6* 

icp 1 — cos 6 i<j> sin 6 i<f> 1+cos 6 

C 2 a/2 2 



(A27) 



XM"9"0) 
X_ 1O (0"£"O) 



V2 



DWe"o) 



cos 6" 



V2 



U) v^ 0"O), 



£)1 ('J."/)" 



and thus, 



(A28) 
(A29) 
(A30) 



(A31) 



and correspondingly 
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|p"10)i = |p"lO) n . (A32) 



For A = 1 we obtain 



X^'WO) = - {(I + cos0cos0') cos(0 - 0') + sin0 sine'} 

-Ucos9 + cos e')sin(0 - 0') (A33) 

X^((j)"e"0) = {- cos9 sin & cos(0 - 0') + sine cose' + i sin e'sin(0 - 0')} (A34) 
v2 

Xl n (0"0"O) = ^{(1 - cose cose') cos(0-0') -sin0 sin^'} 

--(cose- cose') sin(0-0'). (A35) 
2 

Hence, for A = 1 apparently X|, 1 (0"e"O) differs from _D A , 1 (0"e"O), and correspondingly 
|p"ll)jj from |p"ll)j, by a phase factor. Now the difference between Eqs. (A28)-(A30) and 
Eqs. (A33)-(A35) is connected to the value of A. Therefore, the phase factor must depend 
on A and is independent of A'. The latter can be understood as we see that there is no 
A'-dependence in Eqs. (A17) and (A18). In addition the phase factor also depends on the 
set of angles (0, 9, 0', 9'). Thus we write 

Al, A (0"e"O) = e^ A J Di, A (0"e"O), (A36) 

where f2 depends on the set of angles (0, 9, 0', 9'). Noting that £>Q 1 (0"e"O) is real the phase 
Q can be given through its tangential as 

= Jm{X o 1 1 (0"e"O)} = sine'sin(0-0') 

i?e{X o 1 1 (0 // e // O)} — cos e sin 9' cos(0 — 0') + sin 9 cos 9' ' 1 ' 

The Q calculated in Eq. (A37) is also valid for other combinations of A' and A, since Q is 
independent of A' and A. After all these evaluations we summarize that 

4(p')it!s(p)|zSA> = e lQA ^(p")|z5A) (A38) 

DS, A (<f>»9»0) ■ (A39) 

We have restored the spin notation S, since Eqs. (A38) and (A39) are general and hence 
apply to arbitrary spin S, including 5 = 0. 
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FIGURES 



FIG. 1. The spin averaged differential cross section d ^ a de [mb/(MeV • sr)], the analyzing power 
A y and the polarization-transfer coefficient D s i for the (p,n) charge exchange process at projec- 
tile energy = 100 MeV and neutron laboratory scattering angle 6i a b = 13°. The solid line 
represents the 3D calculation, the dashed line the PW calculation with j = 7, J = 31/2. Both 
calculations are based on the Bonn-B potential [17]. 

FIG. 2. Same as Fig. 1, but for projectile energy E[ a b = 197 MeV. The solid line represents 
the 3D calculation, the dashed and dotted lines represent PW calculations with different 2N total 
angular momentum j and total 3N angular momentum J as indicated in the figure. All calculations 
are based on the Bonn-B potential. 



FIG. 3. The spin averaged differential cross section dE a de [mb/(MeV • sr)], the analyzing power 
A y and the polarization-transfer coefficient Du for the (p,n) charge exchange process at projectile 
energy E^ = 197 MeV and neutron laboratory scattering angle 0i a b = 24°. The solid (short 
dashed) line represent the 3D calculations for the first order term based on the AV18 (Bonn-B) 
potential. The long-dashed (dotted) lines stand for the partial wave based, full Faddeev calculations 
based on the AV18 (Bonn-B) potentials. The data are taken from Ref. [17]. 

FIG. 4. The same as in Fig. 3, but for the spin averaged differential cross section d ^ a de 
[mb/(MeV.sr)], the analyzing power A y and the polarization-transfer coefficient D ss at 9^ = 37°. 



FIG. 5. The spin averaged differential cross section dE u dQ [mb/(MeV • sr)] and the analyzing 
power Ay for the (p,n) charge exchange process at projectile energy E\ a i = 100 MeV and neutron 
laboratory scattering angle Q[ a i = 24° The solid (short dashed) line give the nonrelativistic 3D 
calculations for the first order term based on the AV18 (Bonn-B) potential. The long-dashed 
(dotted) lines represent the 3D calculations that include the effects of relativistic kinematics. 



FIG. 6. The spin averaged differential cross section dE a dQ [mb/(MeV • sr)], the analyzing power 
Ay and the polarization-transfer coefficient Du for the (p,n) charge exchange process at projectile 
energy E^ = 197 MeV and neutron laboratory scattering angle 9^ = 24°. The solid (short 
dashed) line give the nonrelativistic 3D calculations for the first order term based on the AV18 
(Bonn-B) potential. The long-dashed (dotted) lines represent the 3D calculations that include the 
effects of relativistic kinematics. The data are taken from Ref. [17]. 

FIG. 7. The same as in Fig. 6, but for the spin averaged differential cross section dE a dd 
[mb/(MeV.sr)], the analyzing power A y and the polarization-transfer coefficient D ss at 6i a b = 37°. 
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FIG. 8. The spin averaged differential cross section ^ [mb/ (MeV • sr)], the analyzing power 
A y and the polarization-transfer coefficient D s i for the (p,n) charge exchange process at projectile 
energy = 346 MeV and neutron laboratory scattering angle 0^ = 22°. The solid (short 

dashed) line give the nonrelativistic 3D calculations for the first order term based on the AV18 
(Bonn-B) potential. The long-dashed (dotted) lines represent the 3D calculations that include the 
effects of relativistic kinematics. The data are taken from Ref. [9]. 

FIG. 9. The spin averaged differential cross section d ^ a de [mb/(MeV • sr)], the analyzing power 
A y and the polarization-transfer coefficient Du for the (p,n) charge exchange process at projectile 
energy E\ a i, = 495 MeV and neutron laboratory scattering angle 9^ = 18°. The solid (short 
dashed) line give the nonrelativistic 3D calculations for the first order term based on the AV18 
(Bonn-B) potential. The long-dashed (dotted) lines represent the 3D calculations that include the 
effects of relativistic kinematics. The data are taken from Ref. [21]. 
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Fig. 4 
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Fig. 8 
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Fig. 9 
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